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Abstract 

Particle packing problems have fascinated people since the dawn of civilization, and continue to 
intrigue mathematicians and scientists. Resurgent interest has been spurred by the recent proof of 
Kepler's conjecture: the face-centered cubic lattice provides the densest packing of equal spheres 
with a packing fraction ip ~ 0.7405 ,l]. Here we report on the densest known packings of congruent 
ellipsoids. The family of new packings are crystal (periodic) arrangements of nearly spherically- 
shaped ellipsoids, and always surpass the densest lattice packing. A remarkable maximum density 
of (/? « 0.7707 is achieved for both prolate and oblate ellipsoids with aspect ratios of \/3 and l/\/3, 
respectively, and each ellipsoid has 14 touching neighbors. Present results do not exclude the 
possibility that even denser crystal packings of ellipsoids could be found, and that a corresponding 
Kepler-like conjecture could be formulated for ellipsoids. 
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BACKGROUND 



Dense packings of nonoverlapping particles have been employed to understand the struc- 
ture of a variety of many-particle systems, including glasses p], crystals heterogeneous 
materials 0], and granular media P]. The famous Kepler conjecture postulates that the dens- 
est packing of spheres in three-dimensional Euclidean space has a packing fraction (density) 
Lp = tt/vTS ~ 0.7405, as realized by stacking variants of the face-centered cubic (FCC) 
lattice packing. It is only recently that this conjecture has been proved ^^O]- Very little is 
known about the most efficient packings of convex congruent particles that do not tile three- 
dimensional space 0, 1^ • The only other known optimal three-dimensional result involves 
infinitely long circular cylindrical particles: the maximal packing density 0max = 7r/\/T2 
is attained by arranging the cylinders in parallel in the triangular lattice arrangement 
Of particular interest are dense packings of congruent ellipsoids (an affine deformation of a 
sphere) with semi-axes a, b and c or, equivalently, with aspect ratios a = b/a and (3 = c/a. 

In two dimensions, it can easily be shown that the densest packing of con grue nt ellipses 
has the same density as the densest packing of circles, if = n / \^T2 ^ 0.9069 [lol This 
maximal density is realized by an affine (linear) transformation of the triangular lattice of 
circles. Such a transformation leaves the density unchanged. In three dimensions attempts at 
increasing the packing density yield some interesting structures, at least for needle-like ellip- 
soids. By inserting very elongated ellipsoids into cylindrical void channels passing through 
the ellipsoidal analogs of the densest ordered sphere packings (an affinely deformed face 
centered cubic or hexagonal close packed lattice), congruent ellipsoid packings have been 
constructed whose density exceeds 0.7405 and approaches 0.7585 in the limit of infinitely 
thin prolate spheroids (ellipsoids of revolution), i.e., when (3=1 and a — > cxo ^, Q. 

However, there appears to be a widespread belief that for nearly spherical ellipsoids the 
highest packing fraction is realized by an affine transformation (stretch by a and jS along 
two perpendicular axes) of the densest sphere packing, preserving the density at 0.7405. 
Mathematicians have often focused on lattice packings, where a single particle is replicated 
periodically on a lattice to obtain a crystal packing. For ellipsoids, a lattice packingis just an 
affine transformation of a sphere packing, and therefore a theorem due to Gauss enables 
us to conclude that the densest lattice ellipsoid packing has ip ~ 0.7405. The next level of 
generality involves nonlattice periodic packings (lattice packings with a multiparticle basis). 
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where a unit cell consisting of several ellipsoids with at least two inequivalent orientations is 
periodically replicated on a lattice to fill Euclidean space. We will refer to these as crystal 
packings. 

We present here a family of crystal packings of ellipsoids that are denser than the best 
lattice packing for a wide range of aspect ratios in the vicinity of the sphere point a = (3 = 1, 
and for certain aspect ratios yields the densest known ellipsoid packings with ^ 0.7707. 



THE DISCOVERY 



We recently developed a molecular dynamics technique for generating dense random pack- 
ings of hard ellipsoids The simulation technique generalizes the Lubachevsky-Stillinger 
(LS) sphere-packing algorithm 1J| to the case of ellipsoids. Initially, small ellipsoids are 
randomly distributed and randomly oriented in a box with periodic boundary conditions 
and without any overlap. The ellipsoids are given velocities and their motion followed as 
they collide elastically and also expand uniformly, while the unit cell deforms to better ac- 
commodate the packing. After some time, a jammed state with a diverging collision rate is 
reached and the density reaches a maximal value. 

Using this technique, we generated nonequilibrium random close packings of ellipsoids 
and observed that for certain aspect ratios (close to a ~ 1.25, (3 ~ 0.8) the random packings 
had a density as high as 0.735, surprisingly close to what we believed was the densest 
ordered packing (stretched FCC lattice) [l5|. This brought into question what the maximal 
density really was for those aspect ratios. Extensive experience with spheres has shown that 
for reasonably large packings, sufficiently slowing down the growth of the density, so that 
the hard-particle system remains close to the equilibrium solid branch of the equation of 
state, leads to packings near the FCC lattice HQ- 

This however requires impractically 
long simulation times for large ellipsoid packings. By running the simulation for very small 
unit cells, from 4 to 16 particles per unit cell, we were able to identify crystal packings 
significantly denser than the FCC lattice, and subsequent analytical calculations suggested 
by the simulation results led us to discover ellipsoid packings with a remarkably high density 
of ~ 0.7707. This result implies a lower bound on the maximal density of any packing of 
congruent ellipsoids, namely, v^max > 0.7707. 
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SUPERDENSE LAMINATED PACKINGS OF ELLIPSOIDS 



We now describe the construction of a family of superdense crystal packings of ellipsoids. 
We start from the FCC lattice, viewed as a laminate of face-centered planar layers of spheres, 
as illustrated in Fig. We similarly construct layers from the ellipsoids by orienting the 
c semiaxis perpendicular to the layer, while orienting the a and h axes along the axes of 
the face-centered square lattice defining the layer, as shown in Fig. In this process, we 
maintain the aspect ratio of the squares of side L of the face-centered square lattice defining 
the layer, i.e., we maintain 

4a 

(1) 



which enables us to rotate the next layer by 7r/2 and fit it exactly in the holes formed by the 
first layer. This two-layer lamination is then continued ad infinitum to fill all space. This 
can be viewed as a family of crystal packings with a unit cell containing two ellipsoids. 

We can calculate the minimal distance h between two successive layers (that preserves 
impenetrability) from the condition that each ellipsoid touches four other ellipsoids in each 
of the layers above and below it. This gives a simple system of equations (two quadratic 
equations and one quartic equation), the solution of which determines the density to be 

167ra/j 

Notice that the axis perpendicular to the layers can be scaled arbitrarily, without changing 
the density. We can therefore just consider spheroids with (3 = 1. The density of this 
crystal packing as a function of the aspect ratio a is shown in Fig. and is higher than 
the density of the FCC sphere packing for a wide range of aspect ratios around the sphere 
point a = /3 = 1, symmetrical with respect to the inversion of a between prolate and oblate 
ellipsoids (we consider the prolate case in the equations in this section). Two sharp maxima 
with density of about 0.770732 are observed when the ellipsoids in the face-centered layers 
touch six rather than four inplane neighbors, as shown in Fig. 12)3, i-e., when L = 2a, which 
gives a = for the prolate and a = l/^/S for the oblate case. These two densest-known 
packings of ellipsoids are illustrated in the insets in Fig. and in these special packings 
each ellipsoid touches exactly 14 neighboring ellipsoids (compare this to 12 for the FCC 
lattice) . 
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THE FUTURE 



There is nothing to suggest that the crystal packing we have presented here is indeed the 
densest for any aspect ratio other than the trivial case of spheres. Many other possibilities 
exist for laminated packings with alternating orientations between layers, and one such 
example is shown in Fig. El More generally, we believe it is important to identify the densest 
periodic packings of ellipsoids with small numbers of ellipsoids per unit cell. This may be 
done using modern global optimization techniques, as has been done for various sphere and 
disk packing problems. However, this is a challenging project due to the complexity of the 
nonlinear impenetrability constraints between ellipsoids. In particular, the case of slightly 
aspherical ellipsoids is very interesting, as the best packing will be a perturbation of the FCC 
lattice with a broken symmetry, and should thus be easier to identify. In Fig. we see that 
the density of our crystal packing increases smoothly as asphericity is introduced, unlike for 
random packings, where a cusp-like increase is observed near a = 1 [15j. Is there a crystal 
packing which leads to a sharp increase in density for slightly aspherical ellipsoids? Further 
multidisciplinary investigations are needed to answer this and related questions. The results 
of such investigations could be used to formulate a Kepler-like conjecture for ellipsoids. 



* S. T., A. D. and F. H. S. were supported in part by the Petroleum Research Fund under Grant 
No. 36967- AC9, and by the National Science Foundation under Grant Nos. DMR-0213706 and 
DMS-0312067. P. M. C. was partially supported by NASA under Grant No. NAG3-1762. 

[1] T. C. Hales. An overview of the kepler conjecture. LANL e-print Archive, 
http : //xxx . lanl . gov/math . MG/981 1071, 1998. 

[2] R. Zallen. The Physics of Amorphous Solids. Wiley, New York, 1983. 

[3] N. W. Ashcroft, N. D. Mermin, and D. Mermin. Solid State Physics. International Thomson 
Publishing, 1976. 

[4] S. Torquato. Random Heterogeneous Materials: Micro structure and Macroscopic Properties. 

Springer- Ver lag. New York, 2002. 
[5] S. F. Edwards. Granular Matter (A. Mehta, editor), chapter The Role of Entropy in the 

Specification of a Powder, pages 121-140. Springer- Verlag, New York, 1994. 



5 



[6] T. Aste and D. Weaire. The Pursuit of Perfect Packing. lOP Publishing, 2000. 

[7] J. H. Conway and N. J. A. Sloane. Sphere Packings, Lattices, and Groups. Springer- Verlag, 

New York, 3rd edition, 1999. 
[8] C. Zong. Sphere Packings. Springer- Verlag, New York, 1999. 

[9] A. Bezdek and W. Kuperberg. Applied geometry and discrete mathematics: DIM ACS Ser. 
Discrete Math. Theoret. Comput. Sci. 4 (P- Gritzmann and B. Sturmfels, editors), chapter 
Packing Euclidean space with congruent cylinders and with congruent ellipsoids, pages 71-80. 
Amer. Math. Soc, Providence, RI, 1991. 
[10] L. Fejes Toth. Regular figures. Pergamon Press, 1964. 

[11] J. Pach and P. K. Agarwal. Combinatorial Geometry. Wiley-Interscience, 1st edition, 1995. 
[12] J. M. Wills. An ellipsoid packing in of unexpected high density. Mathematika, 38:318-320, 
1991. 

[13] A. Donev, S. Torquato, and F. H. Stillinger. An event-driven molecular dynamics algorithm 
for nonspherical particles: Ellipses and ellipsoids. In preparation, 2003. 

[14] B. D. Lubachevsky and F. H. Stillinger. Geometric properties of random disk packings. J. 
Stat. Phys., 60:561-583, 1990. See also Ref.Q- 

[15] A. Donev, I. Cisse, D. Sachs, E. A. Variano, F. H. Stillinger, R. Connelly, S. Torquato, and 
P. M. Chaikin. Jammed ellipsoids beat jammed spheres. Submitted for publication, 2003. 

[16] A. R. Kansal, S. Torquato, and F. H. Stillinger. Diversity of order and densities in jammed 
hard-particle packings. Phys. Rev. E, 66:041109, 2002. 

[17] S. Torquato, T. M. Truskett, and P. G. Debenedetti. Is random close packing of spheres well 
defined? Phys. Rev. Lett, 84:2064-2067, 2000. 

[18] B. D. Lubachevsky, F. H. Stillinger, and E. N. Pinson. Disks vs. spheres: Contrasting prop- 
erties of random packings. J. Stat. Phys., 64:501-525, 1991. Second part of Ref. Q. 



6 



FIGURES 



Figure 1: Part a (top): The face-centered cubic packing of spheres, viewed as a laminate of 
face-centered layers (bottom layer is colored purple and the top layer yellow). Part b (bottom): A 
nonlattice layered packing of ellipsoids based on the FCC packing of spheres, but with a higher 
packing fraction. 

Figure 2: Part a (top): The density of the laminate crystal packing of ellipsoids as a function of 
the aspect ratio a (/? = 1). The point a = 1 corresponding to the FCC lattice sphere packing is 
shown, along with the two sharp maxima in the density for prolate ellipsoids with a = \/3 and 
oblate ellipsoids with a = as illustrated in the insets. Part b (bottom): The layers of the 

densest known packing of ellipsoids, as illustrated in part a. The layers can be viewed as either 
face-centered or triangular. 

Figure 3: Another computer-generated layered packing of prolate ellipsoids (a = \/3) with a zig- 
zag-like variation in the orientations of the ellipsoids in adjacent layers. This packing is slightly 
denser than the FCC lattice {(p = 0.7411), but we did not attempt to find the optimal configuration. 
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Figure 1: Donev et al. 
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